Experimental verification of fidelity decay: 
Prom perturbative to Fermi Golden Rule regime 
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The scattering matrix was measured for a flat microwave cavity witfi classically chaotic dynamics. 
The system can be perturbed by small changes of the geometry. We define the "scattering fidelity" 
in terms of parametric correlation functions of scattering matrix elements. In chaotic systems and 
for weak coupling the scattering fidelity approaches the fidelity of the closed system. Without free 
parameters the experimental results agree with random matrix theory in a wide range of perturbation 
strengths, reaching from the perturbative to the Fermi golden rule regime. 
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The stability of quantum motion has been a topic of 
increasing interest in recent years. In Ref. l], Peres pro- 
posed to consider the time evolution of wave packets gov- 
erned by two shghtly different Hamihonians. Starting 
from the same initial state, their overlap provides a nat- 
ural measure for the stabihty of the quantum evolution. 
As "fidelity" and "quantum Loschmidt echo" , this quan- 
tity has since been investigated extensively (see Ref. Q 
and references therein). Nowadays, it has become a stan- 
dard benchmark for the reliability of quantum informa- 
tion processing Following Ref. [2|, one may define 
fidelity as F{t) = \f{t)\'^ and fidelity amplitude as 



f{t)^{m\uHt) u'{t)\m) 



(1) 



where the unitary operators U'{t) and U{t) describe the 
perturbed and unperturbed time evolution of an arbi- 
trary initial state V'(O)- Depending on the strength of the 
perturbation one can discern three regimes. In the per- 
turbative regime, where time-independent perturbation 
theory can be applied, the decay of the fidelity is Gaus- 
sian. For larger perturbations a cross-over to exponential 
decay is observed, with a decay constant obtained from 
Fermi's golden rule For very strong perturbations 

the decay constant saturates at the classical Lyapunov 
exponent Q. 

Since the first spin-echo experiment by Hahn Q, echo 
experiments have been performed with many different 
quantum and classical wave systems (e.g. Ref. 0, P). 
However, wave functions are usually not accessible to ex- 
periments, and only some reduced information is avail- 
able, such as the nuclear induction averaged over the 
probe in a magnetic resonance experiment jlOl or 
the transmission between two antennas in a microwave 
or ultrasound experiment 0, 0, Q| . 



Here, we report on the experimental measurement 
of fidelity decay in a fiat electromagnetic cavity, using 
the equivalence of Helmholtz and stationary Schrodinger 
equation . Instead of following the evolution of wave 
packets, we measure stationary spectra of scattering ma- 
trix elements, separately, for the perturbed and the un- 
perturbed system. Then, for a given scattering matrix 
element, we compute the Fourier transform of the cross- 
correlation function between the two spectra. After an 
appropriate normalization, this defines the scattering fi- 
delity amplitude. Averaging this quantity over a large 
number of uniformly distributed antennas with small 
transmission yields the standard fidelity amplitude. Yet 
for integrable systems, this may still lead to system spe- 
cific results. For chaotic systems, by contrast, all antenna 
positions become equivalent and measurements for a few 
of them are sufficient. This statement is made more pre- 
cise in jl6j . Its validity is supported by our experimental 
results which agree with the universal prediction of ran- 
dom matrix theory [3, • 

Our microwave experiments are adequately described 
by the statistical scattering theory JJJ , if absorption 
is taken into account The scattering matrix for a 

billiard with two antennas can be written as: 



E — Hcs 

Here Hcs = Hint — {i/2)VV^ is the effective Hamilto- 
nian of the open system and Hmt is the Hamiltonian of 
the closed billiard. The column vectors of V, denoted 
by V^°'\ contain the information on the coupling to the 
antennas at positions Va- For antenna diameters small 
compared to the wavelength, the coefficients Vja are pro- 
portional to il)j{fa), tbe wave functions of the closed sys- 
tem at the antenna positions. 
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Consider the cross-correlation function between an S- 
matrix element of the perturbed and the unperturbed 
system in the time domain: 



(3) 



where S'^^{E) is given by Eq.((21), but with iJint replaced 
by -ff/nf The brackets denote an energy window and/or 
ensemble average. Note that C[S*ij, S'^^]{t) describes a 
kind of echo-dynamics, which is similar to the quantum 
echo defined in Eq. , but decays even without any per- 
turbation. We therefore use the autocorrelation function 
for a heuristic normalization and define the scattering 
fidelity amplitude as 



fab{t) = 



m 



(4) 



Reflection and transmission measurements have been 
performed in a flat microwave cavity, with top and bot- 
tom plate parallel to each other . The cavity is quasi- 
two-dimensional for frequencies v < c/{2h), where h is 
the height of the billiard. The billiard is shown in the in- 
sert of Fig.n ft consists of a rectangular cavity of length 
L = 438 mm, width B = 200 mm and height h = 8 mm, 
a quarter-circle insert of radius i?i = 70 mm, and a half- 
circle insert of radius i?2 — 60 mm placed on the lower 
side. The position of the latter was changed in steps 
of 20 mm to generate an ensemble average over 15 dif- 
ferent systems. Additional elements were inserted into 
the billiard to suppress bouncing-ball resonances: two 
half-circle inserts with radius i?3 = 30 mm, and a slope 
on the upper boundary. The perturbation of the sys- 
tem was achieved by varying the length L in steps of 
I = n ■ 0.2 mm, with n — 1-10. The change of area and 
surface due to the shift of the billiard wall, was taken into 
account by unfolding the spectra to a mean level distance 
of one. The frequency window of the Fourier transforms 
was 1 GHz wide, and a Welch filter was applied. In this 
range the antenna coupling and the wall absorption are 
approximately constant. 

We compare the experimental results with the random 
matrix prediction 



f{t) = exp 



l-JJUr')dT'dr 



(5) 



This expression is obtained by exponentiation of the lin- 
ear response result, thus incorporating the known behav- 
ior in both, the perturbative and the Fermi golden rule 
regime |l7l |. Here, A is the perturbation strength, and 
&2(i) is the two-point form factor for the Gaussian or- 
thogonal ensemble. We use dimensionless units, where 
the Heisenberg time tn = h/A is equal to one and A 
denotes the mean level spacing. 

Figure shows on a logarithmic scale the cross- 
correlation function C[S*i, S[i]{t) given in Eq. to- 
gether with the autocorrelation function C[S*i, Sii]{t). 
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FIG. 1: (Color online): Logarithmic plot of the correlation 
function C[Sii, S'n] for u — 5-6 GHz, I = 1mm and A — 
0.047. The experimental results for the auto correlation are 
shown in black, while the correlation of perturbed and unper- 
turbed system are shown in gray / orange. The smooth solid 
curve corresponds to the theoretical auto-correlation function, 
and the dashed curve to the product of auto-correlation func- 
tion and fidelity amplitude. The insert shows the billiard 
geometry used. Moveable parts are marked with an arrow. 
See text for dimensions. 
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FIG. 2: (Color online): Logarithmic plot of the scattering 
fidelity amplitude fii{t) for v = 5-6 GHz, I = 1mm and A — 
0.047. The smooth curve shows the linear-response result for 
the fidelity amplitude f{t), where the perturbation strength 
A was obtained from the variance of the level velocities. 



The latter agrees with the corresponding theoretical au- 
tocorrelation function, calculated in 19]. The parame- 
ters for the wall absorption and the coupling of the anten- 
nas have been obtained according to the same reference. 
With increasing time, C[Sli, S[i]{t) deviates more and 
more from the autocorrelation function. This deviation 
contains the essential information on echo dynamics. 
For the frequency range shown in Fig. ^ the pertur- 
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bation strength A was determined directly from the mea- 
sured spectra via the variance of the level velocities. The 
dashed curve in Fig.^is a product of the theoretical auto- 
correlation function and the fidelity amplitude (Eq. 
of the closed system. The experimental result for the 
cross-correlation function (Eq. ©) agrees perfectly over 
six orders of magnitude with the linear-response expres- 
sion without any free parameter. This justifies our defini- 
tion of the scattering fidelity amplitude fab{t) in Eq. 10} . 

FigureElshows fab{t), computed from the experimental 
cross- and autocorrelation functions according to Eq. Q). 
There are two advantages for using the experimental and 
not the theoretical autocorrelation function: The com- 
puted quantity does not depend on theoretical assump- 
tions, and the influence of non-generic features, visible in 
the correlation functions, is reduced. 

We now study the dependence of the fidelity decay on 
the perturbation strength. In our experiment A varies 
from A = 0.01 for n = 1 and v = 3-4 GHz up to A = 0.5 
for n = 10 and = 17-18 GHz. Figure 13 shows the 
scattering fidelity amplitude for three different frequency 
windows. Here, A has been fitted to the experimental 
curves, as its determination from the level dynamics is 
time consuming and for strong perturbations not always 
feasible. To improve statistics, experimental results for 
/iij /22 and fi2 have been superimposed. 

For the random matrix model one expects a transi- 
tion from linear to quadratic decay near the Heisenberg 
time. In the perturbative regime, the linear term in the 
exponential is still close to one and we observe Gaus- 
sian decay of the fidelity amplitude, as seen in Fig.l^Ia). 
With increasing perturbation strength the linear term 
becomes more pronounced, leading to the Fermi golden 
rule regime The (exponentiated) linear-response 

formula ij^l agrees very well with experiment throughout 
the range. Recently, an exact solution for the random 
matrix model proposed in [T^ has been obtained using 
supersymmetry techniques 21]. This result is shown as 
dashed lines in Fig. For the accessible perturbation 
strengths the experiment does not allow to distinguish 
between the linear-response and the exact result. 

In billiard systems, the parameter variation is not due 
to a change of the Hamiltonian, but of the boundary 
condition. It was shown in chapter 5 of Ref. that both 
situations are equivalent. For the case of a parameter 
variation in the billiard due to a shift of a straight wall 
the matrix element of the equivalent perturbation reads 



(Hi 



d'ipn{x,y) difjm{x,y) 



dx 



dx 



dy, (6) 



x=0 



where I is the shift of the wall (in x-dircction) and B is 
the length of the shifted wall. The perturbation strength 
according to the random matrix model in [T^ is given by 
the variance of the off-diagonal matrix elements: 



(7) 
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FIG. 3: Logarithmic plot of the fidelity amplitude f(t) for 
three different perturbation strengths. Linear response result 
(smooth solid gray line), exact theoretical result (dark dashed 
line), and experimental result (solid dark line). The perturba- 
tion parameter A has been fitted to each experimental curve. 
The parameters were v — 3-4 GHz, / — 0.4mm, A = 0.01 
for (a), u = 13-14 GHz, / = 0.6mm, A = 0.13 for (b), and 
V = 16-17 GHz, I = 0.8mm, A = 0.21 for (c). 



Using Berry's conjecture of the superposition of plane 
waves we can derive for large wave numbers k: 



2L 
3^ 



3i2 



3^ 



3;2 



(8) 



In Ref. 23] , the same expression has been obtained using 
periodic orbit theory and the ergodicity assumption. 
Figure^shows the experimental perturbation strength 
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FIG. 4: (a) Perturbation strength as a function of the 
shift I of the billiard wall for the frequency windows v = 'i- 
4 (stars), 9-10 (diamonds) and 16-17 GHz (triangles). The 
slope of the straight lines is 2. (b) as a function of the 
frequency v for I = 0.2 (stars), 0.6 (diamonds) and 2.0mm 
(triangles). The slope of the straight lines is 3. 



as a function of the shift I of the billiard wall for 
three different frequency regimes (a), and as a function 
of the frequency, for three different shifts (b) . We observe 



excellent agreement with the scaling A 



predicted 



in Eq. 0. This can be used to average all experimental 
data as a function of the scaled variable 47r^A^C(i), thus 
reducing the fluctuations almost completely . In spite 
of the correct scaling, the experimental prefactor is about 
three times smaller than predicted. In cases where we 
determined the variance of level velocities directly from 
the measured spectra, we found the same discrepancy. 
The deviation is caused by the fact that we are far from 
the semiclassical limit, for which Eq. (jSj) was derived. 
Additional numerical studies for the Sinai billiard by H. 
Schanz 24] substantiate this explanation. 

We have shown that the fidelity amplitude fab{t) of 
scattering matrix elements Sab is an easily accessible 
quantity. In our experiment it approximates the ordi- 
nary fidelity amplitude f{t). The relevant parameters 
are well controlled. This enabled us to verify the theo- 
retical results for fidelity decay covering the range from 



the perturbative to the Fermi golden rule regime without 
any free parameter. For strong perturbations the exact 
solutions of random matrix theory predict a revival of 
the fidelity amplitude at the Heisenberg time, whereas 
semiclassics predicts an exponential decay with the Lya- 
punov exponent up to the Ehrenfest time . It remains 
an open question, whether these features can be verified 
in a future microwave experiment. 
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